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Asymptotics of an optimal compliance-location 

problem 

G. Buttazzo* F. S ant a mb r o g i o tin d N. Varchon^ 


Abstract: We consider the problem of placing a Dirichlet region made by n 
small balls of given radius in a given domain subject to a force / in order to mini¬ 
mize the compliance of the configuration. Then we let n tend to infinity and look 
for the T—limit of suitably scaled functionals, in order to get informations on the 
asymptotical distribution of the centres of the balls. This problem is both linked to 
optimal location and shape optimization problems. 
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1 Introduction 

The study of asymptotical problems in the optimal location of an increasing amount 
of resources has been developed intensively, even in recent times, mainly by using an 
approach based on T—convergence. In |4] the so-called location problem (choosing 
a set £ composed by n points in a domain O in order to minimize the average 
distance of the points of hi from £) is investigated as n —> +oo, finding a T—limit 
of a suitable sequence of functionals on the space of probability measures on fl 
In jjlBj the same analysis is performed for the so-called irrigation problem, where 
points are replaced by connected dimensional sets of finite length, and the constraint 
jj£ < n by 7d 1 (£) < l. Both the problems are linked to the Monge-Kantorovich 
optimal transport theory. However, these asymptotical problems are not completely 
understood since explicit minimizing sequences are not in general known, apart some 
simple cases, usually in dimension two. For instance, for the location problem it is 
known that placing the points on a regular triangular grid, so that each one is in the 
middle of a cell shaped like a regular hexagon, gives an asymptotically minimizing 
sequence (see H2j or M for stronger results). 

On the other hand, many researches have been carried out on shape-optimization 
problems involving PDEs, i.e. optimizing the shape of a domain where to solve a 
PDE (usually of elliptic type with prescribed boundary conditions), in order to min¬ 
imize the value of an objective functional depending on the solution of the PDE. 
There is a wide literature on shape optimization problems, both from a theoretical 
and numerical point of view. We refer for instance to the books CQ. [2.i :5J. 
ra, where the reader can find various approaches and a lot of examples and details. 
It is well known (see for instance [BJ and [7]) that for general cost functionals the 
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existence of an optimal domain may fail and a relaxation procedure, involving the 
use of capacitary measures, is needed. It is also known that for cost functionals ful¬ 
filling particular additional monotonicity assumptions, a simple volume constraint 
is sufficient to imply the existence of a minimizer (see for instance L S|). Moreover, 
these problems have shown to have many industrial applications in engineering and 
mechanics (bridges, light structures supporting loads...) and this is the reason 
for their widespread study and especially for numerical computations. One of the 
simplest shape optimization problem, which is also one of the most important in 
applications, is compliance minimization. It consists in finding a domain O (usually 
under a volume contraint) which minimizes the compliance value f (} fun, dC d where 
Uq is the solution of the elliptic equation —A u = f with Dirichlet boundary condi¬ 
tions on dfl. Such a problem satisfies the monotonicity assumptions needed to have 
the existence of minimizers. 

What we consider in this paper (see Section □ is a compliance minimization 
problem where the unknown domain where to solve the PDE with Dirichlet boundary 
conditions is the complement of a finite union of balls whose number and radius are 
assigned. In fact, for given n G N and r n > 0, we look at the problem of choosing 
n balls {B(xi,r n ))i = \ t .., in in order to minimize the compliance of f2 \ UT=i B(xi,r n ) 
(Q C and / G L 2 (Q) are fixed). Obviously, the problem is meaningful only under 
a condition like r n < crT x ^ d , otherwise the total volume of the balls may be sufficient 
to cover all 0, thus obtaining a vanishing compliance. 

The aim of the present paper is to look at the optimal location of those balls 
when n —> +oo. Using small balls is the most natural way to approximate a point- 
location problem, which would not make any sense since points have zero capacity in 
dimension two or more, and so they do not affect the solution of elliptic equations. 
In this paper we just consider the case r n = an -1 ^, for a fixed small parameter 
a, which corresponds to a volume constraint. Moreover, in Section Owe deal with 
the one dimensional case, getting it as a particular case of the d— dimensional, but 
highlighting the case of points instead of balls as well (which corresponds to taking 
a = 0). 

The domain has to be considered as an elastic membrane and / stands for 
the forces acting on it (for instance the different loads it has to carry). Our goal 
is to reinforce or support the membrane at some points (by mechanical devices, 
by some kind of glue...), choosing where to locate the support points, so that it 
bends as less as possible, minimizing the work of the forces, given by f fu dx. The 
important issue here is that we let the number of support points increase, reducing 
correspondingly the effect at each point, and we look at the asymptotic density of 
support points. 

As it happens for the location problem, placing n balls with given radii in Q 
reduces to a finite dimensional variational problem as a consequence of the severe 
geometric constraints imposed to the admissible domains. Our asymptotic analysis 
is based on some techniques developed in Jj and mainly in ra and so this paper, 
even if dealing with a shape optimization subject, stands apart from other shape 
optimization papers both for the setting of the problem and for its development. 

Since we solve elliptic PDEs in open domains which are obtained by removing 
small holes from a fixed one, this subject is in connection with the problem of homog¬ 
enization in perforated domain (see for instance PI), Yet, there are some important 
differences. First, we have an optimization problem on the holes instead of taking 
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them as given. Consequently, no periodic assumptions on the holes is supposed, 
even if periodic structures are used many times in the proof of the T—limit result. 
Finally, homogenization on perforated domains is interesting when the size of the 
holes decreases quicker than n -1 / d , otherwise, under Dirichlet boundary conditions, 
the solutions trivially tend to 0. This is the reason for introducing a scaling factor 
which enlarge the values of the functionals and of the solutions when n increases. 

An interesting question we do not consider in the present paper is how to deal 
with the case r n = a n n~ x l d with a n —> 0, since this could be considered as a better 
approximation of the case of points. Moreover, there are further reasons to study 
in the future this new case. First, the case a n —> 0 is more linked with the theory 
of perforated domains developed in Pi. where the radius of the ball tends to zero 
faster than n~ 1 ^ d . Then, this case seems to require better mathematical techniques 
involving the behaviour of the Dirichlet energy when radii tend to 0, since it is not 
simply possible to build a recovery sequence by homogenizing a given configuration 
(since this would let the total volume of the balls unchanged, while here it has to 
tend to zero). Finally, the result we get here ('Theorem 12.21) involves a function 9 
we are not able to compute explicitly, while the case a n —> 0 seems to require just 
its asymptotic expansion near 0. 


2 Locating balls to optimize compliance 

For any open set D C M d , a > 0 and neNwe define: 


A(a, 7i)(Q) 


sen 


n 

D fl B(xi,r ) for x* G D r , r = arT 1 ^ 1 

Z— 1 




where stands for the r—neighbourhood of f2. Given D C and / G L 2 (fl), for 
any compact set E C 11 with positive measure let us define the function Mys.o as 
the solution in the weak sense of the problem 


—An = / in D \ S 
u = 0 in S U dfl, 


which means precisely u G Hq {Q \ E) and 


[ Wu-V^dC l = [ f(j)dC d for any 0 G (2.1) 

Jo. Jn 

Notice that / > 0 implies M/,s,o A 0, by the maximum principle. For f > 0, we 
define the compliance functional over subsets of a given domain D as 

F(E,/,D)= [ fu f ^ndC d = [ \ Vu /iS , n | 2 dC d . 

Jn Jn 

Remark 1. The requirement / > 0 seems to be mostly a technical assumption to 
simplify the proofs (it allows us to deduce some pointwise inequality by maximum 
principle on the solutions of Dirichlet problems): physically speaking it means con¬ 
sidering only forces which have the same direction on the whole D (for instance usual 
gravity). 
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What we want to do now is considering such a compliance functional on the set 
•4.(a,n)(f2) (notice that imposing E G A(cx,n)(D) implies a volume constraint given 
by |E| < w<ioi d and a geometrical constraint, i.e. compelling E to be composed by an 
assigned number of identical balls). This is our n —th compliance minimization prob¬ 
lem. The following existence result holds and can be proven by standard methods, 
due to the very severe geometry the elements of A(a,n)(Q) are constrained. 

Theorem 2.1. For any n G N, if Q is any bounded open subset ofM. d and f > 0 
belongs to L 2 (Q), the problem 

min{F(E,/,fi) |E G A(a,n)(n)} (2.2) 


admits a solution. 


Then we would like to let n tend to infinity and look at the asymptotics of the 
problem, mainly at the distribution of the centres of the balls. Let us associate to 
every E G A(a, n)(Q) a probability measure on fi, given by /iv. = ra -1 Y^=i 
where p : W l —> is a fixed projection of the whole space to fi. The role of the 

projection p is simply to handle the case where the centre of the ball B(xi,r ) lies 
outside fi. Such a measure is an atomic measure uniformly distributed on the centres 
(or on their projections). Then we define a functional F n : Tift) —> [0, +oo] by 

F( , = f n 2/d F(E,/, Q) if p, — /i£, E G A(a, n)(0); 

^ ] Too otherwise. 


The coefficient n 2 ^ d is a factor which is needed in order to avoid the functionals to 
degenerate to the trivial limit functional which vanishes everywhere. Anyway, such 
a coefficient does not affect the choice of the minimizers. 

We will give a T—convergence result for the sequence ( F n ) n , when endowing the 
space V{Fl) with the weak* topology of probability measures. To introduce the limit 
functional F we need to define the quantity: 

9{a ) := inf {lirninf n 2/d F(E n , 1 ,I d ) |E n G A(a,n)(I d ) } , (2.3) 

where I d = (0, l) d is the unit cube in W l . This quantity will play the role of the 
constant appearing both in ra (as 9 n}P ) and in [3] (as Cf)- However, its dependence 
on a will be essential and, unfortunately, in general not explicit. It is easy to see 
that 9 is a decreasing function on M + , which vanishes after some point. In fact 
if a > \fdj 2, it is possible to use n balls of radius an~ l ^ d to build a set E G 
*4.(a, n)(I d ) covering the whole cube F l , thus getting a vanishing solution Uf ^,i d = 0 
and F(E, /, I d ) = 0. Let us call t\ the first vanishing point, i.e. 

t 1 := inf (f G K 1 9(t) = 0} < 

We denote by 9~ and 9 + the lower and upper semicontinuous envelopes of 9 , respec¬ 
tively. They are given by 


9 ( a ) 
9 + (a) 


sup { 9(/3 ) |(3 > a} 

inf {9((3) \/3 < a } . 
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It is easy to check that the following formula holds: 

0~(a) = inf jliminf n 2 //d F(E n , 1 , I d ) |S n G A(a n , n)(I d ), a n —» a j . (2.4) 

We may now define the candidate limit functional F by setting, for /i G 7^(0) 

F(p)= [ (2.5) 

J [Aa 

where /x a denotes the density of the absolutely continuous part of g with respect 
to the Lebesgue measure. It is evident from m that the whole behaviour of the 
function 6 affects the minimization problem for F. It is more convenient to introduce 
the function g a defined by g a {x) = x^ 2 ^ d 9(ax 1 / d ), which is a decreasing function 
whose semicontinuous envelops gf and g~ are obtained by < 7 +(a;) = x~ 2 / d Q+ [ax l ^ d ) 
and g“(x) = x^ 2 ^ d 9~(ax 1 ^ d ). So we have 

Jn 

Here is a sketch of the behaviour of g a , according to what already highlighted 
and to what proven in Section 0 on the function 6. In particular g a is a convex 
function. 



Figure 1: Qualitative behaviour of the function g a 

The result we will prove is the following: 

Theorem 2.2. Given any bounded open set C a non-negative function f G 
L 2 (Q) and a > 0, the sequence of functional ( F n ) n previously defined V—converges 
towards F with respect to the weak* topology on V{VL). 

The consequences of such a T—convergence result, by means of the general theory 
(see mi), are the following: 

• for any sequence (E n ) n of optimal sets for the minimization problem (12.211 it 
holds, up to subsequences, gz n — 1 p where // is a minimizer of F\ 

• should F have a unique minimizer, we would get full convergence of the whole 
sequence n^ n to the unique minimizer /i. 

• the sequence of the values inf {F(E, /, O) |E G A(a, n)(D,) } is asymptotical to 
n~ 2 / d inf {F(/j) \/i G V(kl) }. 
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It turns out that it is very important to investigate about the minimizers of the 
limit functional F, which is much related to the behaviour of g. From Section 0 we 
know that g a = g~ and that g a is a convex and strictly decreasing function (up to 
the point t x where 6 = 0). If we define t Q = ti/a 1 it holds g a (x) = 0 for any x > t a . 
We will restrict our analysis to the case where a is sufficiently small, so that it is 
not possible to cover the whole by n balls of radius cm~ 1//rf . We summarize in the 
following statement what we know on the minimizers of F. 

Theorem 2.3. If a is such that a < any minimizer g for F is an absolutely 

continuous probability measure with density g a which satisfies 

Pafx) G {-dga)- 1 (j^j) , 

for a suitable constant c > 0. In particular it holds g a <t a . 

Proof. Being g a strictly decreasing up to t a it is straightforward that optimality 
implies absolute continuity: otherwise, just remove the singular part from g and use 
the same mass on the absolutely continuous part, enlarging its density and strictly 
decreasing the value of F. Moreover, by using Lagrange multipliers or performing 
simple variations to g. it is easy to get the existence of a constant c > 0 such that 

f 2 (x) {-dg a (fj, a (x))) 3 c, (2.6) 

where dg is the, possibly multivalued, subdifferential of g. Then we get 

g a {x) G (-c ^)- 1 > 

where, for a multifunction G, we use the notation G^ 1 (t) = {z G [0, +oo] : t G G(z)}. 
We sketch a possible behaviour of —dg a and (— dg a )~ l in Figure EH 



Figure 2: The sub differential g' a of g a and its inverse 


To get g a < t a (which is suggested by Figure [21 as well) it is sufficient to notice that 
otherwise in (El one should have c = 0. This would imply g a > t a a.e. and this is 
not possible by the assumption on a, since g has to be a probability measure. □ 

A consequence of g a < t a is the fact that strict convexity of g a on (0, t a ) is 
sufficient to ensure uniqueness of the minimizers of F. Section 0 will also provide 
some qualitative property of 9, in order to guess the behaviour of g a . 

6 








3 The T —convergence result 

We will prove Theorem 12.21 in several steps, the most important two corresponding 
to the T — lim inf and T — lim sup inequalities. 

Proposition 3.1. Under the same hypotheses of Theorem denoting by F~ 
the functional T — lim inf n F n , it holds F~(p) > F(p) for any p G 'P(fl). This 
means that, for any sequence (£ n ) n such that weakly--k converges to p and E n G 
*4.(a, n)(fi), it holds lim inf n n 2 ^ d J fu n dC d > F(p), where u n stands for Uf^, n ,a- 

Proof. First of all, let us fix e > 0 and, in analogy to what performed in Ha. define 
the set G s>n as follows: let us assume C [—a, a] d for a G N, then we set 

Ge,n = (J B(y,r), r = an ~ 1/d , k = [{en ) 1,d J . 

y^k~ 1 Z d n[— a,a] d 

Now we define E' n = U G Etn and we set u' n = Uf t w n ,n- Since u n > u' n , it is 
sufficient to estimate from below the integrals n 2//d J fu' n dC d . The utility of the 
new sequence (u' n ) n lies in the fact that (n 2 ^ d u' n ) n is L 2 — bounded. In fact we have 
0 < u' n < Uf t G £:n ,n and, by Lemma HOI it holds 

||n/,G e .n,dU 2 (n) < C(a, e, f)n 2/d . 

This implies that (7i 2 ^ d u' n ) n is bounded in L 2 and so, up to a subsequence, we have 
rUl d u! n w. Since in this case we have lim inf „ n 2 ^ 1 j fu n dC d > j fwdC d , it is 
sufficient to estimate w from below. To do this, we first estimate the average of w 
on a cube Q centred at point x G 12. It holds 

[ w dC d — limn 2/Vi [ u' n dC d . 

Jq n Jq 

We use 

u' n > Uf^Q = Uf( x ) t v n ,Q + Uf-f( x ),V n ,Q > Uf( X ),V„,Q - ^ Q, 

where the first inequality comes from the fact that we add Dirichlet boundary con¬ 
ditions on Q and the last by maximum principle. Notice that this is the key point 
in the proof where we strongly use f > 0 (in fact the other pointwise inequalities 
could be replaced by global integral estimates valid for general /, but here this is 
not possible, since we have to estimate integrals performed on Q instead of on fl). 

We will estimate separately the two terms. Let us start from the easiest, i.e. 
the latter. It holds < u \f-f(x)\,G £ , n ,Q- By applying Lemma O to the 

domain Q we get 

\\u\f-f(x)\,E' n ,Q\\LPQ < n' 2,d C(a,e)\\f - f(x)\\ L 2 ( Q y 
Then we estimate, by Holder inequality, 

n2/d [ u \f-fU)\,T.’ n ,QdC d < n 2r/ \Q\ 1/2 \\u lf _ f{x)]x , ntQ \\ L 2 {Q 

Jq 

< C{a,£)\Q\ 1/2 \\f - f( x )\\ L 2 {Q) . 
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(3.1) 



We now evaluate the other term. First we define the number 


k(n, Q) = j} (ji : B(xi, an l / d ) n Q ^ 0 j U | j : B(yj , an x / d ) D Q ^ 0 j j . 

Now notice that m/( x ),s^,q = f{^)ui,E' n ,Q- Let us denote, for simplicity, the functions 
Ui t z' t Q by v n . By a change of variables, if A is the side of the cube Q and we define 
v n ,\(x) = X~ 2 v n (Xx), it holds v U) \ = u lx -±Y,' n j d - We notice that 

Moreover, it holds k(n, Q ) —»■ +oo, since 


k(n,Q) > j : B(y j ,an~ 1 / d ) nQ ^ 0j « en\Q\. (3.2) 

We may also estimate the ratio between k(n, Q ) and n, by using (13.211 and the 
fact that $ < i : B(xi, an~ 1 / d ) n Q ^ 0 > = //. n {Qs n ), where 8 n = an~ x ! d —> 0 and Q$ 


denotes the 8 —neighbourhood of Q. From /i n —*■ // we have limsup n i^n(Qs n ) < l^{Q), 
so that 

limsup ^ < /x(Q) + e|< 5 |- 


n 


(3.3) 


Now, by using the equality in d2HD, it is not difficult to derive that 

liminf k(n, Q) 2 ^ d [ v n dC d = lirninf k(n, Q) 2 ^ d X d+2 [ v n \dC d 

n Jq n Ji d 


> »"r (jWg)+£|oi) 1/d ). (3,4) 


So we get 


liminf n 2//d [ v n dC d = lim ( ° ] liminf k(n,Q) 2 ^ d [ v n dC d 

n J Q n \k(n, Q) J 


> A««-(H0, W ) +e |Q|)^( 


<2 

1 


^ rl “(?§ +£ 


1 /dN 


I <31 


/i(Q) + e|< 3 | 

f Jgl A 

\d{Q) e \Q\) 


2 /d 


2 /d 


This implies, recalling also CD. 


I QQ wdC“ > -C( a ,e)\Q\-^\\f-f(x)\\ L , m 
Jq 

+ Jgl 

V \Q\ J ) \h(Q)+ £ IQI 


2 /d 


Now we let Q shrink towards x, thus obtaining, using lower semicontinuity, for a.e. 
x G 


w(x)>f(x)6 (a(fj, a (x) +e) 1,d 


Ha(x) + £ 


2 /d 















where we have used the fact that a.e. point x is a Lebesgue point for / G L 2 to 
get the first term vanish, and that for any measure y and a.e. x 6 0, it holds 
\Qpy(Q) —> /i a (ic) when the cube Q shrinks at its centre x. So we get 


1 im inf n 2/d / fu n dC d > 


> 


fwdC d 

r 


(a {n a {x) + e) 1/d ^j dC d , 


In (da(x) + e) 

and our original aim is achieved when we let e —> 0, i.e., still using that 9~ is l.s.c., 

P 


liminfn 2/d / fu n dC d > 


In Ta\ x 


y/d 


9 {ay a (x) 1 l d ') □ 


Lemma 3.2. The following facts hold. 

1. For any 0 < £q < 1 there exists a constant C = C(e o) such that 


v G H\l d ), | {n = 0} | > e 0 | I d 


' I d 


v 2 dC d < C / |Vn| 2 d£ a 

Jld 


2. If we replace I d by a cube Q whose side is A the same is true with the constant 
A 2 C instead of C. 

3. As a consequence, for any e > 0, any n G N, any domain and f G L 2 (Q) 
with F > 0, the function Uffj € n ,n, where G e ^ n is defined as in Theorem, \3 . 11 
satisfies ||M/,G e ,„,n||L2 ( n ) < n~ 2/d C(a, £)\\f\\ L 2 {n) . 

Proof. The first assertion comes from a well-known variant of Poincare inequality 
and can be proven by contradiction. The second one is just a scaling of the first. 
To prove the last, let us consider a family of cubes Qj, each centred at a point 
yj G /c _1 7L d n [—a, a] d , whose side is k~ l . Let us extend the function Uf t G e , n ,n to the 
set O' = U. Qj , where the union is over the cubes touching by defining a function 
v which is identical to Uf t a e , n ,n on ^ and 0 outside (we recall that we have Dirichlet 
boundary conditions on so that such an extension belongs to H l [yt')). Notice that 
v vanishes in a whole ball of radius an~ l ! d in each of these cubes, so that the ratio 
between the volume of such a ball and the volume of the cube depends only on a 
and £ and not on n. By applying the second part of the statement of this lemma, 
we get 

[ v 2 dC d <C{a,e)n~ 2/d I \Vv\ 2 dC d , 


Qj 


Qj 


and, by summing over j, we get 

[ v 2 dC d <C{a,£)n~ 2/d [ |Vu| 2 dC d . 


in' 


in 1 


Since v vanishes outside we may write 

[ dC d < C(a,E)n- 2 ' d I |V« /l0 ,.. 


'dC° 


fuf'Ge.n&dC 1 < C(a, e)n 2/d ||u/,G e ,„,n||L 2 (n)||/|U 2 (n). 


The thesis easily follows by dividing by | |u/,G e , n ,n| U 2 (n)- 


□ 
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To get also the opposite inequality, i.e. the T — limsup inequality, we need this 
crucial lemma. 

Lemma 3.3. Given E 0 G A(a 0 , n 0 )(I d ), a domain 0 C M d and f G L 2 (12), we 
consider the sequence of sets 

E fc = |J (y + k-'Eo) HO. 

y£ k~ 1 Z d 

It happens E fc G A(ao, n(k, 12))(12), where n(fc, 12) m \I}\k d n 0 . Then we consider the 
sequence ( Uk)k, given by 

Uk = k 2 Uf t syn- 

If we assume dl d C E 0 , then it holds Uk —*■ c(E 0 )/ where the weak convergence is in 
the L 2 sense and c(E 0 ) is a constant given by f Jd u 1So jd dC d . 

Proof. First, we notice that the sequence (uk)k is bounded in L 2 (12). This may be 
proven in a way very similar to that of Lemma 13.21 but here it is even simpler. In 
fact we may extend as before Uk to the union of the cubes of the kind y + k~ 1 I d 
which intersect 12, with y G /c _ 1 Z d , by giving it the value 0 outside 12. Then we apply 
standard Poincare inequality to each cube Q y (Uk vanishes on the boundary of Q y 
since dl d C So), and here the Poincare constant is k~ 2 C , where C depends only on 
the dimension d. Then, by putting all the inequalities together, restricting to 12 and 
integrating by parts as in Lemma 1X21 we get ||xtfc[|z, 2 (f 2 ) < C , ||/||l 2 (q). Let us now 
consider an arbitrary weakly convergent subsequence (not relabeled) and its limit 
Wf ) s 0 ,o- It is easy to see that the pointwise value of this limit function depends only 
on the local behaviour of /. In fact, the key assumption dl d C E 0 produces small 
cubes around each point x G 12 which do not affect each other. So, if / = JT /j/^ is 
piecewise constant (the pieces Ai being disjoint open sets, for instance), it happens 
that for large k the value of Uk at £ G Ai depends only on /). It turns out that, 
for a piecewise constant function /, it holds Wf t s 0 ,o = fwi,z 0 ,n- It is indeed clear 
that in this case (/ = 1 ), since we are simply homogenizing the function w 1 So /d, the 
limit of the whole sequence {uk)k exists, does not depend on the global geometry 
of 12 , but it is a constant and it is the same constant as if there was I d instead of 
12. Then the constant is easy to be computed and is the constant c(E 0 ) appearing 
in the statement. It remains now just to show that the equality Wf t z 0j n = fc( E 0 ) is 
true for any L 2 function /. The convergence of the whole sequence will then follow 
easily by uniqueness of the limit of subsequences. To get the result for a generic /, 
just take a sequence (/„)„, of piecewise constant functions approaching it in L 2 and 
notice that 

k 2 Uf,xi<,n = k 2 u fn ^k : n + /c 2 M/_/ ni sfc,o- 

The first term here weakly converges to f n c( Eo) as k —> +oo, while the second is 
bounded in the L 2 norm by C\\f — f n \\L 2 {Q.)- This means that any weak limit of 
subsequences of (uk)k must be close in the L 2 norm to f n c( E 0 ), i.e. 

\\ w f, So, O - fnC{ E 0 ) | |z,2(fl) <C\\f— f n 11 L 2 (n), 
which implies, letting n —> +oo, s0 ,o = fc( E 0 ). □ 

Now we want to build efficient sets Eo satisfying the key assumption of our 
previous Lemma, that is dl d C Eo (we will call those sets for which such an inclusion 
boundary-covering sets). 
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Lemma 3.4. For any a > 0 and any £ > 0 there exists no € N such that for any 
n > tiq we find a" < a and a set Eg A(a" ,n)(I d ) which is boundary-covering, with 
n 2 / d f Id u i^jddC d < (1 + e)9 + (a). 

Proof. Let us fix <5 > 0 and a' < a such that O(a') < (1 + 5)9 + (a). By definition of 
9(a'), we may find e A(a', ni)(I d ) such that 


n 


2/d 

1 



Ui jS i id dC d < (1 + S)9(a') 


and, moreover, the number n\ may be chosen as large as we want. Now we enlarge 
the set Si to get a new set E 2 which is boundary-covering: we add to some 
m balls of radius r = a'nf l ^ d (the same radius of the balls composing Zb). In 
order to cover dl d the number of balls we need does not exceed C fr d ~ x , so we have 
m < C(oi')n\ _1//d . It is possible to choose rii so that 


m < Sn i and a 


^ni + m'j 1//d 


= a 


< a. 


This is useful, since 


S 2 G A , n\ + m \ ( I d ) = A(a" ,n 2 )(I d ), 

where we set n 2 = n± + m. Moreover 

n^ d [ u lT , 2 iddC d < (—^ n 2 / d [ /d dC d < (1 + S) 2+2 ^ d 9 + (a). 

Ji d ’ \ n \ J Ji d 

Now, if we are given a large number n, we just need to homogenize the set E 2 . 
By homogenization of order k of a set S C I d into a domain A we mean the set 
A D |J yG k ~i Z d y + k-'S. Here we take the homogenization of order k of S 2 into I d , 
which is a set 

E e A{a\ k d n 2 )(I d ) c A(a", n)(/ d ), 

where we choose k such that k d n 2 < n < (k + l) d n 2 . For this set E it holds 
( k d n 2 ) 2 / d f Jd Ui ^ jddC d = ruj d [ Jd Uip 2 jdd£ d , so that 

n 2,d £ Ul ^jddC d < (^T^J ‘ (1 + S) 2+2/d 9 + (a). 

If n > n 2 ( 1 + <5 _1 ) d , then it holds k > 5 _1 and 1 + 1/k < 1 + S, so that we get 

n 2/d [ Ul s !ddC d < (1 + 5) 2+A/d 9 + (a). 

Jld 

It is now sufficient to choose S sufficiently small so that (1 + S) 2+4 ^ d < 1 + e and 
then set n 0 = n 2 ( 1 + <5 _1 ) d . □ 


11 



We are now ready to start with the T — lirnsup main part. We will start from 
a very particular class of measures. Let us call piecewise constant those probability 
measures p G V(Ut) which are of the form 

p = p ■ C d , with p G L l (Ul), / pclC d — 1, p > 0, 

Jn 

for a piecewise constant function p = YiT=o Pi^hi the pieces hi, being disjoint Lips- 
chitz open subsets with the possible exception of h2 0 = \ IJ”=i To simplify the 

notation, let us also define the functional F, which is the same as F with the only 
difference that we replace 6~ by 9 + : 

F(P)= [ f 9+{a d ] dC“ = f fgMdC d . 

Jn pi in 

Proposition 3.5. Under the same hypotheses of Theorem \2.A it holds 
F + (p) < F(p), where F + = F — liiri sup F n , 

n 

for any a piecewise constant measure p G V{UL). This means that, for any such 
a measure p and any e > 0, there exists a sequence of sets (E n )„ such that p^ n 
weakly-* converges to p, £ n G A(a,n)(kl) and moreover it holds 

limsup n 2/d [ fu f ^ n ^dC d < (1 + e) [ f 2 ° ^ K lC d . 
n Jn Jn P ' 

Proof. First of all, let us consider the numbers ap 1 ^, which will appear as arguments 
of the function 6 + . We know by applying Lem ma id. 41 to all of them that there exists 
a common number no and some sets £* G A(a" p]/ d , no)(/ d ), which are all boundary¬ 
covering and such that 


r h) d [ u i s i iddC d < (1 + e)6 + (api). 

Jv> 

Now, if we are given some numbers k(n,i), we define the sets E^ by homogenizing 
into the set £, of order k(n,i ) i.e. 


^ = U y + k(n,i ) x £\ 

y&An.i 

where A Uii = {|/G k(n,i)~ 1 Z d \(y + k(n, i) _1 £*) D ^ 0}. We define as well 

K= (J y + KnA)- 1 ^, 

y£An,i 

i.e. without intersecting the balls with flj. 

Then we choose £ n = |J- TJ n U £ n where E n is a union of balls of radius a"n~ l ^ d 
covering the union of the boundaries dkli inside the interior of hi. The number of 
balls we use in E n is approximatively \kli\k(n, i) d no in each zone hlj, plus Cn 1 ~ l / d for 
the set E n , where C depends on a" and the total perimeter of the partition In 
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each zone 0* the radius of the balls is given by a"pl^riQ 1 ^/k(n, i). This means that 
the sequence £ n we are building is admissible (i.e. £ n e .A(a,n)(fl) and /is„ —*■ //) 
if we have 


l/d 


, ^ 1 -< an 1//rf for n large enough and for i — 0 ,..., m; 

n 0 k(n,i ) 


\Qi\k(n, i) d n 0 + CA 1 l ^ d < n and is asymptotic to n; 
k(n, i) d n 0 


i =0 


n 


Pi for i = 0,..., m. 


The hrst conditions involves radii, the second one the number of balls and the third 
one is related to weak convergence. All these conditions are satisfied if we set 


k(n, i ) = 


(<i “ 


l/d 


We want now to estimate the values of the functionals on this sets E n . We have 
used the set £ n covering the internal boundaries of the sets in order to get a local 
behaviour in which different zones Q,: are independent on each other. The quantity 
we want to estimate is 

P m r, m p 

n2/d / f u f,x n ,n dC d = fn 2/d u f ,x n ,n dC d < / fn 2/d u f ^ Ui dC d . 


In the last inequality we have used the fact that £ n allows adding Dirichlet boundary 
conditions on each fh and T, l n D H l n . The disintegration of the integral here per¬ 
formed allows us applying Lemma l3.3l on each Q,. Notice first that n ks k(n, i) d p~ 1 no. 
Then, by Lemma HOI we know that k(n, i) 2 Uf^ n ,nt —*■ c(£*)/, where the convergence 
is weak in L 2 . By our choice of £* we have c(£ l ) < (1 + e)n 0 2,/d 6 l+ (ap i 1 ^), so we get 

lim n 2/d [ fu fA ^.dC d < (1 + £)p~ 2/d 9 + (ap] /d ) f f 2 dC d , 

71 JCli JQ.i 


and, summing up, 


hmn 2/d 


f u f, 




dC d < 



2 9 + (ap 1/d ) 

1 f?/d 


d £ d , 


that is our thesis. □ 

We must extend our result to non piecewise constant measures and replace 6 + 
by 9~. To simplify the notation, we will use the function g a defined in Section [21 

Proposition 3.6. For any p e V{Fl) it holds 

F+(p)<F(p)= [ f 2 g- a (p a )dC d . 

Jn 
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Proof. Given p — p ■ C d , with p > c > 0, we may choose a sequence p n = p n ■ C d 
where p n —> p a.e. and p n are piecewise constant with p n > c. This may be done 
by approximating p in L 1 first by regular functions, then by functions which are 
constant on cubes, for instance. So, by lower semicontinuity of F + (see [IT]), we get 

F + (p) < liminfF(/i n ) < [ f 2 gf (p) dC d , 

11 Jn 

where we used the inequality gf(p n ) < < 7 +(c) to have dominated convergence of 
9a(Pn) t° their a.e. limit, which is estimated by upper semicontinuity by gf(p). So 
we have extended our inequality to any absolutely continuous measure with positive 
density bounded away from 0. Now take p = p ■ C d without the assumption p > c 
and take p c = (p + c)1 p <m + (p — £ c )l p >m, where e c is chosen so that f n p c dC d = 1. 
It is clear that e c —* 0 when c —> 0 and that p c = Pc' £ d converges to p. So we have 

F + {p) < lim inf F(p n ) = lim f f 2 gf (p c ) clC d + lirn f f 2 gf(p c )dC d 

c ^° c ^° Jp<M c ^° Jp>M 

= [ f 2 g~ (p) dC d + [ f 2 gf (p) dC d 

J p<M J p>M 

< [ f 2 g-(p)d£ d + gf(M) [ f 2 d£ d . 

Jq J p>m 

The number M being arbitrary, we get easily F + (p ) < f Q f 2 g~(p) dC d by letting 
M tend to +oo. We need now only to get the result for measures of the form 
p = p ■ C d + v with a singular part v. Such a measure may be approximated by a 
sequence of measures p n = plA n + vl where the v— part carries the mass of the 
singular part v and 1a„ —>> 1q in L 1 . We may choose v so that it is v > M with 
arbitrary large M, and in particular we may have gf(v) = 0 since the function 6 
vanishes for large values of the argument. So we get 

F + (p) < = lim inf [ f 2 g~(p)dC d < [ f 2 g~(p)dC d . 

71 n J A n Jn 


With this the proof is over. 


□ 


4 Some properties of the function 0 

In this section we study some properties of the function 6. We already know that it 
is decreasing and that it vanishes from a certain point on. Our first goal is to show 
the equality 6 = 9~, so that we can get rid of the l.s.c. envelopes. 

Proposition 4.1. For any a > 0 it holds 9(a) = 9~(a). 

Proof. From its definition: 

9(a) = inf |lim inf n 2 /d F(S n , 1 , Q) |E n e A(a, n)(Q) j 
= lim inf n 2 ^ d min {F(E, 1, I d ) |E 6 A(a,n)(I d ) \ . 

n L 1 J 


14 


It is then clear, by T—convergence, that such a lim inf is in fact a limit equal to the 
minimum of the limit problem in I d with f — 1. So 


6(a) = min / g a (p) dC a 


' i d 


p > 0, p G L\l d ), / pdC d = 1 

JI* 


since the minimum will be certainly achieved on an absolutely continuous measure 
(as already mentioned in Theorem 12.31 as a consequence of g~ being decreasing). 
Moreover, it is clear that p = 1 achieves the minimum: in fact, taken an optimal p, 
we may homogenize it so that we get a sequence p n weakly converging to the constant 
1. Each of these measures will be a minimizer since the functional only takes into 
account the measure of the level sets of p n , which are the same as in p. So, by lower 
semicontinuity, p = 1 will be a minimizer. We deduce 6(a) = g~ (1) = 6~(a), and 
this yields that 6 is lower semicontinuous. □ 

Now we have g a = g&'■ our next step will be the following: 

Proposition 4.2. For any a > 0, the function g~ is convex. 

Proof. This is a consequence of the lower semicontinuity of F, which is implied 
by the fact that it is a T—limit (see mi ). This is quite standard (by necessary 
conditions on functionals of this kind in order to be l.s.c., see for instance [3]), but 
not immediate since we are restricted to probability measures. Let us prove it, for 
the sake of completeness. In fact, given two values x, y > 0 and t G [0,1] it is possible 
to build a probability measure p on a domain R = (0, a) x I d ~ l with density x in 
(0, at) x / d_1 and y in (at, a) x F l ~ l (a has to be chosen so that a(tx+ (1 — t)y) = 1). 
By homogenizing p on the Erst coordinate we get a sequence p n in which the value 
of the functional F is constant and equal to F(p) (we take / = 1). Then we 
notice p n —^ p^, where p^ is the measure with density tx + (1 — t)y. By lower 
semicontinuity we get 


ag (tx + (1 - t)y) < atg (x) + a(l-t)g (y), 
which gives the required convexity. 

Corollary 4.3. The function 6 is locally Lipschitz continuous on (0, +oo). 


□ 


Proof. This is a consequence of the equality 6(x) = gi(x d )x 2 , where both the factors 
in the right hand side product are locally Lipschitz functions of the variable x. □ 

The next thing to do with 6 is to prove some estimate. In particular it is necessary 
to get estimates from below, so that 6 is not identically 0, otherwise our limit 
functional F would be trivial. Anyway, we start from an estimate of 6 from above: 

Proposition 4.4. There exists a constant C = C(d) such that, for any a > 0, it 
holds 6(a) < Ca 2 ~ d , for d > 3, or 6(a) < Clog((v / 2o:)~ 1 ), for d = 2. 

Proof. It is sufficient to consider a particular sequence (or even a subsequence) of 
sets £ n G A(a, n)(I d ) and then to compute the lim inf in the definition of 6(a). 
We choose to consider just the numbers n of the form n = k d , and to build, for 
each k G N, a set which is composed by n = k d balls of radius a/k, with their 
centres placed at the middle points of the k d cubes of side 1/k of a regular lattice 
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partitioning the cube I d . First, we notice that it holds id < v n , where v n is the 
solution to the problem 


-Av n = 1 
V n = 0 

Fk l ’ n = ^ 


in I d \ S n 
in E n , 
on dl d . 


The inequality between the solutions of Dirichlet and Neumann problems comes from 
maximum principle. By scaling arguments, it is clear that the energy n 2 A f Jd v n d,C d 
equals the energy f Jd V\ dC d . If we set r'o = \fd/ 2, which is the radius of the smallest 
ball containing the cube I d and centred at its same centre, it holds V\ < w, where 
w is the solution of the Neumann problem on such a ball: 


{ -Aw = 1 in B(x 0 , r 0 ) \ B(x 0 , a) 
w = 0 in B(x o, a), 

= 0 on dB(xo, r 0 ), 

where Xq is the centre of the cube I d . This solution may be explicitly computed, it 
is radially symmetric and is given by 

w(x) = k(d)(a 2 ~ d — r 2 ~ d ) — \(r 2 — a 2 ) with r — \x — x 0 \, if d > 3; 

2d 

w(x) = k log ( — ) — -(r 2 — a 2 ) with r = \x — x 0 |, if d — 2. 

\aJ 4 

It turns out that 

^-w > 0 on dl d , 
on 

which allows using maximum principle to get V\ < w. So it is sufficient to compute 
the integral of w on J d , which can be estimated by k(d)a 2 ~ d in the case d > 3, and 
by /dog(a _1 ) + ATogr 0 = A; log((v / 2ct) _1 ) in the case d — 2. □ 


A similar estimate from the other side may be obtained, as we show in our next 
proposition. Here the techniques we use are quite different and much more related 
to shape optimization and PDEs. 


Proposition 4.5. For a <t\ the following estimates hold: 

„1 -d 


-o' {a) > 

0(a) > 
0(a) > 


a 


- 2 - 
dud d 

a 2 ~ d 

d(d — 2 )u>d 
log(a -1 ' 1 


2t r 


— C if d > 3, 
C if d = 2, 


(4.1) 

(4.2) 

(4.3) 


for a suitable constant C = Cfd), where u>d denotes the measure of the M d —unit ball. 


Proof. For a fixed n G N and fixed points (xi ) i= G I d , let us consider the sets 
= IXiS^.on- 1 /-) G A(a, n)(I d ) and Q a = I d \ E a . The following estimate 
holds, by Holder inequality: 



I 9 I 


1 dH 


d -1 


> 



a ,i d dH 


d -1 


io*r 


(4.4) 
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where the last equality comes from integrating by parts J n —A u 1Sa jd dC d . Notice 
that it holds as well 


H d 1 {d£L a ) < ndu d (an l ^ d ) d 1 = du d a d 1 n 1//d . 
Then, we recall that, by using shape derivative (0), it holds 


d 


d 


—F(E a ,l,I d ) = n- 1 ' d I \—u 1A:I <\ 2 dH i - 1 , 


da 


I an a 


(4.5) 


(4.6) 


since we perturbate by a vector held which is normal to the boundary of the 
balls and proportional to n~ 1//d . So we have, by putting together m, EHD and 


<nni) 


^n 2 / d F(E Q ,l,/ d )> 

dot du)d 


a d 1 |Q a | 2 a d 1 (1 — 2c v d a d ) 


> 


du d 


where in the last inequality we have used |il a | 2 > (1 — to d a d ) 2 > 1 — 2uj d a d . So, by 
integrating over an arbitrary interval (aq,a!o) C (0,£i), we have 


n~ 


2 / d F(Z ai ,l,I d )>n 2 / d F(Z ao ,l,I d ) + 


pao 


du d 


a 1 “(1 — 2u d a d ) da. (4.7) 


'ai 


We now compute explicitly the right hand side of dUZD in the case d > 3, getting 


2 -d „,2 -d 


n 


2/d F(i: ai , 1 , I d ) > n 2 / d F(Y> ao) 1 , I d ) + -± 


a i —a 


«q — a 2 


d{d — 2)u d d 
Then, passing to the inf over ( Xi)i and to the lirninf over n, we get 


9(a i) > 9(a 0 ) + -± 


2 -d 2-d 

Oi- i Q(n 


0^0 Oi.-^ 


(4.8) 


d{d — 2)u d d 
This gives the estimate on 9 we were looking for, with 

C(d) = t\/d + t\~ d /[d{d - 2)u! d ]. 

Anyway, from (EE1). we can also infer the estimate on 9': we let an —> a 0 and divide 
by a 0 — aq. On those points «o where 9 is differentiable (i.e. almost everywhere, 
since 9 is locally Lipschitz) it holds exactly 

-«'M > ~ 2 V- 

dujd d 

In the case d = 2 it is sufficient to compute again the integral in EH), getting 
nF( E«, 1. I d ) > „F(E O0 ,1. I d ) + M^I -logK- 1 ) _ oj-af _ 

Z7T Z 

Then the conclusion follows in the same way. Here C = log(t7 1 )/27r + t\/2. □ 

We summarize now all the results on 9 we have got in this section. 
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Theorem 4.6. The function 9 :]0, +oo[—> [0, +oo[ is a nonincreasing function, lo¬ 
cally Lipschitz on ]0, +oo[, such that x i—► 9(x)x~ 2 / d is a convex function. Moreover 
9 is not identically 0 but it vanishes from a certain point on, i.e. 9(x) =0 for any 
x >t\ with t\ < y/d/ 2. Finally it holds 

C\x 2 ~ d - C < 9{x) < C 2 x 2 ~ d if d> 3 
Ci | log a;| — C < 9{x) < C 2 | log rr| if d = 2. 

We terminate this section by stressing the interest in finding explicit minimizing 
sequences for the case / = 1 and D = I d , since this would give the value of 9. 
In analogy to what happens in the location problem (see [H; and we may 

conjecture that minimizing sequences are given by placing the centres of the balls 
on some kind of regular grids and, moreover, as far as d = 2, hexagonal regions 
with balls in the middle seem to be good candidates. This would lead to a better 
knowledge of 9 , thus letting us get better information on the minimizer /j for the 
limit functional F by Theorem 12.hi As a weaker conjecture, we may think that the 
function 9 , for which we have proven upper and lower estimates in term of x 2 ~ d or 
| log(x)|, is such that there also exist the limits 

lim ^ for d > 3 or lim ^ for d = 2. 

x->0+ X 2 d a;^0+ | log(x)| 

In fact what we have already proven is just that the ratios above are bounded. For 
instance, in the case d > 3, should 6{x) actually behave like x 2 ~ d for x near 0, we 
could derive, for small a , a behaviour like 

ha ~ cf. 

5 The one dimensional case 

In the case of dimension 1 we are able to compute the function 9 explicitly. 

Everything is, in fact, simpler in dimension 1, since the balls we remove are 
intervals which disconnect the domain of the differential equation (an ODE in this 
case), and so we can compute explicitly the solution. We have already pointed out 
that in dimension 1 the compliance optimization problem is well-posed also for finite 
unions of points (i.e. the case a = 0), and not only for small intervals. So we will 
take into account also the value of 0(0). 

Theorem 5.1. For any a > 0, if d — 1, it holds 

6(a) = (® (1 - 2 “ )3 !/ “-f 

[0 if a > 2 ■ 

Proof. Given an interval J whose length is l > 0, the solution of the problem 

f —u" = 1 in J 
) u = 0 on dJ 

is given by u(x) = x{l — x)/2 (if we suppose J = (0,/)), and so its integral on J 
is Z 3 /12. When we are given a G [0,1], we take into consideration disjoint union 
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of n intervals, for a total length of 2a (in this case it is necessary to have disjoint 
intervals if we want to minimize compliance). So the energy of the configuration is 


if 


under the condition l t 


%— 1 


i= 1 


1 - 2a. 


By convexity of l t—> Z 3 , the minimum of such a quantity is achieved by U = n 3 (1 — 
2a), and so it holds 


0(a) = lirninf n 2 n- - = —(1 — 2a) 3 . 

v ; n 12n 3 12 v ' 

The case a > 1/2 follows trivially from 6( 1/2) = 0 and 9 being non negative and 
decreasing. □ 

As a consequence, the limit problem in dimension 1 is completely known. Up 
to inverting the resulting function g' a it is possible to find explicitly the minimizer, 
thanks to the expression for p, given by Theorem 12.31 

We can also consider the sequence of minimization problems with n points instead 
of n balls: we have to consider 


KM 


n 2 F(S, /, I) if p = /i S and jjE < n 
Too otherwise. 


Theorem 5.2. Given J = [a, b] and f G L 2 (J), the sequence of functionals (F n ) n 
overV(J) T— converges, with respect to weak* convergence onV(J), to the functional 
F given by 


Hh) = 


P 


12 


/ // 2 


dC 1 . 


Moreover, F has a unique minimizer n op t , given by 


hopt = cf 2/i ■ C l , where c= ( / / 2/3 dC 1 


-l 


Proof. This statement follows by slight modifications of what proven in last section. 
The only point where a > 0 was actually used was Lemma fd.21 Anyway, in the one 
dimensional case, the inequality holds if we replace | {f = 0} | > £o|^ d | by {v = 0} ^ 
0. This allows performing again the proof, by following the same steps. To find 
explicitly the minimizer it is sufficient to use Lagrange multipliers, getting that 
/ 2 //i 3 is constant. □ 
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